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Fluctuation-dissipation (FD) relation of the three-dimensional Heisenberg spin glass with weak 
random anisotropy is studied by off-equilibrium Monte Carlo simulation. Numerically determined 
FD ratio exhibits a "one-step-like" behavior, the effective temperature of the spin-glass state being 
about twice the spin-glass transition temperature, T^s — 2Tg, irrespective of the bath temperature. 
The results are discussed in conjunction with the recent experiment by Herisson and Ocio, and with 
the chirality scenario of spin-glass transition. 
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Off-equilibrium dynamics of spin glass (SG) has at- 
tracted much recent interest [1]. In equilibrium, there 
holds a relation between the response and the correlation 
known as the fluctuation-dissipation theorem (FDT). In 
off-equilibrium, relaxation of physical quantities of SG 
depends on its previous history. Most typically, it de- 
pends not only on the observation time t but also on the 
waiting time tw, exhibits aging. While in the short- 
time quasi-equilibrium regime << t << tu, (to is a 
microscopic time scale) , relaxation is still stationary and 
FDT holds, it becomes non-stationary and FDT is bro- 
ken in the long-time aging regime t >> In particular, 
the breaking pattern of FDT, and its possible relation to 
the static quantity, has been the subject of recent active 
studies. 

A quantity playing a central role here is the so-called 
fluctuation-dissipation ratio (FDR) X, which may be de- 
fined by the relation. 
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where R(ti,t2) is a response function measured at time 
t2 to an impulse field applied at time ti, C{ti,t2) is a 
two-time correlation function in zero field at times ti 
and t2, and T is the bath temperature. One may re- 
gard T/X = Tcff as an effective temperature. In the case 
FDT holds, one has X = 1 and T^s = T. 

Via the study of certain mean- field models, Cuglian- 
dolo and Kurchan showed that, in the limit of infinite 
time ii,t2 oo, the FDR X depended on the times 
ti and t2 only through the correlation function C{ti,t2), 
i.e., X{ti,t2) = X{C{ti,t2)) [2]. It was further suggested 
that X{C) could be related to an appropriate static quan- 
tity, i.e., X{C) is equivalent to the a::((7)-function describ- 
ing the replica-symmetry breaking (RSB) pattern in the 
Parisi's scheme, which is related to the overlap distribu- 
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tion function via P{q) 
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[2]. This conjecture was 



supported both by numerical simulation [3] and analytic 
work [4] (see, however, also ref. [5]). 

Experimental studies on the FDR of SG have been 
hampered for a long time by the difhculty in perform- 
ing high-precision measurements of correlations (noise 
measurements), although some interesting preliminary 



results were reported [6]. Recently, however, a remark- 
able experiment by Herisson and Ocio for an insulating 
Heisenberg SG CdCr1.7Ino.3S4 has eventually opened a 
door to experimental access to the FDR of SG [7]. Some 
comparison was already made between these expeirmen- 
tal results and the numerical results obtained by off- 
equilibrium simulation on certain SG models. 

Meanwhile, it should be remembered that all numer- 
ical results of the FDR so far available are limited to 
the Lsing-like SG models, which represent an extremely 
anisotropic limit [3,5,8-11]. By contrast, the SG material 
studied in ref. [7] was more or less Heisenberg-like. Re- 
cent experimental and numerical studies suggested that 
both equilibrium and off-equilibrium properties might 
differ significantly between the Ising and Heisneberg SGs 
[12-17]. Under such circumstances, in order to make a 
direct link between theory and experiment, it is clearly 
desirable to study the Heisenberg-like SG where magnetic 
anisotropy is weak. The purpose of the present Letter is 
to perform extensive off-equilibrium simulations on such 
a weakly anisotropic Heisenberg-like SG model, and nu- 
merically investigate its FDR. Interestingly, the calcu- 
lated FDR of the three-dimensional (3D) Heisenberg SG 
exhibits a feature of the so-called one-step RSB. 

I consider the weakly anisotropic classical Heisenberg 
model whose Hamiltonian is given by 



n 



(2) 



where Si ={Sf , Sf , S^) is a three-component unit vec- 
tor at the i-th site, and the sum runs over all nearest- 
neighbor pairs of a 3D simple cubic lattice with N ~ 
spins. .Jij is the isotropic random exchange coupling tak- 
ing the value J or —J with equal probability (±J or 
binary distribution), D^^ (Mj ^ = 2;,y, z) is the random 
exchange anisotropy which is assumed to be symmetric 



Z?j ■ ) and uniformly distributed between the in- 



terval [—D,D], and /i is a magnetic field applied along 
z-direction. Aimed at clarifying the properties in the 
limit of weak anisotropy, I assume here a rather weak 
anisotropy of D/J = 0.01. The SG transition temper- 
ature of this model estimated from equilibrium Monte 
Carlo (MC) simulation is Tg/J ~ 0.21 [18]. 
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Dynamical Monte Carlo simulation is performed 
according to the standard single spin-flip heat-bath 
method. At time zero, I quench the system from an infi- 
nite temperature to a working temperautre T below Tg, 
and wait in zero field during the time t^. In calculating 
the response, a weak DC field of intensity h is applied, 
and subsequent growth of the magnetization per spin, 
[< ruzit + tw) >], is recorded at time t + t^t, (< • • ■ > 
and [■ • ■] are thermal average and average over quenched 
randomness). In calculating correlation, I measure the 
two-time correlation function of magnetization at times 
tw axidt+tyj, C[t]tw) = [< 'mz{tw)-'mz{t+tw) >], always 
staying in zero field in this case. 

The integrated FD relation eq.(2) yields the zero- field 
cooled (ZFC) susceptibihty x=[< iTiz{tw + t) >]/h as, 

X{t;t^) = ^f X{C')dC', (3) 

where the factor i takes care of the number of compo- 
nents of Heisenberg spins. 

The lattice size mainly studied is L = 32 with periodic 
boundary conditions. I generate total of 10® Monte Carlo 
steps per spin (MCS) in each run, with = 2.5 x 10^ or 
10^ MCS. Sample average is taken over 40-160 indepen- 
dent bond realizations. 

The calculated spin autocorrelation functions are 
shown in Fig. 1 as a function of observation time scaled by 
the waiting time, t/t^, for several temperatures below Tg. 
The decay of the autocorrelation at lower temperatures 
occurs in two steps. Near the end of the quasi-equilibrium 
regime, the decay curves tend to level off exhibiting a 
plateau-like structure, which is expected to correspond to 
the static Edwards- Anderson order parameter ^ea [12]. 
Furthermore, at lower temperatures T/ J — 0.10 and 0.05 
the two curves for different t.^j cross around t/t.^ ~ 1, 
and a,t t > i^,, the data for larger lie above the one 
for smaller (superaging) . Such a superaging behavior 
means that the aging is apparently more enhanced than 
the one expected from the naive t/tu,-scaling. The ob- 
served superaging of spin correlations is similar to the one 
previously observed in the Heisenberg SG with the Gaus- 
sian coupling [12]. At higher temperatures T/J — 0.15 
and 0.20, by contrast, the curve for longer lies slightly 
below the one for shorter tyj (subaging), or nearly comes 
on top of it for longer t (near full aging) . To examine the 
robustness of the observed aging behavior against the 
system size and the waiting time, I show in the inset the 
autocorrelation for the smaller size L — for the tem- 
peratures T/J = 0.15 and 0.05, with = 5 x 10'^ and 
2.5 X lO'*. The superaging/subaging tendency observed 
for L = 32 can also be seen here. 

Hence, whether the spin correlation exhibits superag- 
ing or subaging apparently depends on the bath tempera- 
ture, superaging at lower temperatures and subaging (or 
near full aging) at higher temperatures. In the present 
case, the borderline temperature is around T/J ~ 0.15, 
corresponding to ~ 0.7Tg. In the recent noise experiment 



by Herisson and Ocio, subaging behavior is observed at 
T ~ O.STg [7]. It might be interesting to see experimen- 
tally whether the superaging behavior as observed here 
ever arises at lower temperatures. In real experiment, 
however, care needs to be taken to an inevitable cooling- 
time effect which might give some bias toward apparent 
subaging [19,20]. 
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FIG. 1. Decay of the spin autocorrelation function of the 
weakly anisotropic 3D Heisenberg spin glass with _D/ J = 0.01, 
plotted versus \og-^^Q(t / t-u,) at several temperature below Tg. 
The lattice size is L = 32, averaged over 80 {T/J — 0.20), 
160 (T/J = 0.15) and 40 (T/J = 0.10) samples. Inset: The 
autocorrelation decay for the smaller system of size L = 16 
averaged over 40 samples. The waiting times are ~ 5x 10'^ 
(blue) and 2.5 x 10* (pink) here. 

In the lower panel of Fig. 2, I show the growth of the 
ZFC susceptibility Tx(t; t^) as a function of t/t^ at sev- 
eral temperatures. I applied a very week field of intensity 
h/J = 0.01. A word of caution is in order here. Gen- 
erally, applied fields should be sufficiently weak so that 
one is in the linear response regime. In the previous sim- 
ulations on the Ising SG, the field intensity of h/J — 0.1 
turned out to be weak enough [3,5,9]. In the case of 
Heisenberg SG, however, still weaker field appears to be 
necessary to guarantee the linearity of the response. This 
is illustrated in the upper panel of Fig. 2 where I show the 
t-dependence of the ZFC susceptibility at T/J = 0.15 
for three different field intensities, h/J = 0.05,0.02 and 
0.01. Some deviation is discernible between these curves 
at longer i, suggesting that the requirement of the lin- 
earity is rather severe in the Heisenberg case. In the 
following, I show the response data for the weakest field 
studied, h/J = 0.01, where the linearity seems almost 
satisfied in the investigated time range. 

Indeed, such a strong nonlinearity of the response ap- 
pears to be even more enhanced in the isotropic limit 
D = 0. 1 also made a preliminary simulation on the fully 
isotropic Heisenberg SG with D = at a temperature 
T/J = 0.15 to study its response, with varying the field 
intensity H/J = 0.05, 0.02, 0.01 and 0.005. While, for a 
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field of moderate intensity H/ J = 0.05, the response x{t) 
of the isotropic system Z? = is almost common with the 
one of the weakly anisotropic system D/J = 0.01, on de- 
creasing the field intesnity, of the isotropic system 
constantly gets larger at longer observation time beyond 
the corresponding value for D/J — 0.01, without show- 
ing any tendency of saturation up to H/J — 0.005. This 
strongly suggests that the field scale associated with the 
linear response is closely correlated with the strength of 
the mangetic anisotropy D, although its detailed proper- 
ties are yet to be clarified. 
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FIG. 2. Growth of the ZFC susceptibility x(t\U?) of the 
weakly anisotropic 3D Heisenberg SG with D/J = 0.01. The 
lattice size is L = 32, averaged over 80 (T j J = 0.20), 160 
(T/J = 0.15) and 40 (T/J = 0.10) samples. Upper panel: 
x(i;iiD) plotted versus the observation time t for several ap- 
plied field intensities. The temperature is T/J = 0.15, with 
the waiting time — 10^ MCS. Lower panel: Tx{t;tm) plot- 
ted versus logjQ(t/tiu) at several temperatures below Tg, with 
U = 2.5 X 10'* and lO'^ MCS. 

As can be seen from the lower panel of Fig. 2, the curve 
for longer tends to lie slightly above the one for shorter 
tw (subaging) , whereas for i > > the two curves for dif- 
ferent merge asymptoticaly (full aging). The superag- 
ing behavior, which has been seen in the correlation, is 
not clearly visible here in the reponse. Experimentally, 
response of SG usually exhibits a weak subaging [1,7] 
or near full aging [20], which seems consistent with the 
present result. 

In Fig. 3(a), I show a parametric plot Tx(t;tw) versus 
C{t; tin) for several temperatures. The straight line in the 
figure is the FDT fine, Tx^ il-C)/3. In the short-time 



quasi-equilibrium regime, the data lie on the FDT line, 
and deviate from the line in the long-time aging regime. 
In the aging regime, the data for longer tend to lie 
above the data for shorter t^, slightly, but systemati- 
cally. Such a systematic drift of the data with varying 
is a common aspect observed both in experiment [7] and 
in simulations on the Ising SG [3,5,8-11]. 

In Fig.3(b), the data for longer are replotted in the 
form of x{t',tw) versus C{t;tu])- One sees that the data 
in the aging reime at various temperatures lie on top of 
each other, indicating that x is temperature independent 
in the aging regime, apparently satisfying Parisi- Toulouse 
hypothesis [3,10]. 
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3. Tx versus C plot (a), and x versus C plot (b), of 



the weakly anisotropic 3D Heisenberg SG with D/J = 0.01 
at several temperatures below Tg. The applied field intensity 
is h/J — 0.01. The lattice size is L = 32, averaged over 80 
(T/J = 0.20), 160 (T/J = 0.15) and 40 (T/J = 0.10, 0.05) 
samples. In (b), tw is lO'' MCS. The broken lines represent 
the FDT lines. The solid line in (b) is the straight-line fit to 
the data in the aging regime, its slope being equal to 0.8. 

Remarkable in Fig. 3 is the strikingly linearity of the 
plot in the aging regime, although some curvature is dis- 
cernible in the part where the data begin to deviate from 
the FDT line. In particular, the data in the aging regime 
appear to come to the C — axis with a nonzero tangent. 
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Such a behavior seems distinct with the behavior of the 
Ising SG where the data exhibit some more curvature, 
and appear to come to the C = axis almost horizon- 
tally, or at least be consistent with a vanishing tangent 
at C = [3,9,10]. By making a linear fit of the data in 
the aging regime in Fig.3(b), I estimate the slope to be 
about 0.8, which corresponds to the effective temperature 
Teff ~ 0.42 J ~ 2Tg, irrespective of the bath temperature 
T(< Tg ~ 0.21J). 

Such a one-step-like behavior of the FDR, though fa- 
miliar in structural glasses, may sound surprising in SG. 
However, this observation finds a natural explanation in 
the recently developed chirality scenario of SG transi- 
tion, which ascribes the true order parameter of the SG 
order to the chirality [12,21]. In this scenario, chirality 
exhibits a one-step-like RSB in the ordered state (chiral- 
glass state) [13], which is manifested in the SG order in 
real Heisenbcrg-like SG via the weak random magnetic 
anisotropy. Thus, the chirality theory predicts the oc- 
currence of a one-step-likc RSB in the SG ordered state 
of real Heisenberg-like SG, which, in its off-equilibrium 
dynamics, might well be reflected in a one-step-like be- 
havior of the FDR. In particular, a nonzero tangent is 
expected at C = for the FDR, which corresponds to 
the delta-function peak in the static overalp distribution 
function P{q) at g = [13]. The present numerical ob- 
servation seems fully consistent with such an expectation 
from the chirality scenario. 

Now, I wish to compare the present numerical result 
of the FDR with the recent experiments. The one-step- 
like linear behavior observed in the present simulation 
seems at least not inconsistent with the experimental 
data of refs. [6,7], although the data of ref. [7] appear 
to exhibit some curvature showing a tendency to level 
off toward C = 0. It would be interesting to experimen- 
tally explore further deep into the aging regime (smaller 
C regime) , to see whether a finite tangent persists or not 
at C ^ 0. Meanwhile, the experimentally determined 
effective temperautre of the SG state of CdCr1.7Ino.3S4 
was Teff ~ mK ~ l.QTg {Tg ~ 16.2K) at the measuring 
temperature T = I'i.ZK [7], which is in good agreement 
with the present estimate of Toff given above. 

It should be noticed that CdCr1.7Ino.3S4, though ba- 
sically a Heisenberg-like magnet, might possess a bit 
stronger anisotropy than the very weak anisotropy as- 
sumed here {D / J = 0.01). Such Ising-like tendency of 
CdCr1.7Ino.3S4 might be a possible cause of the curva- 
ture observed in the x-C plot of ref. [7]. It would then 
be highly interesting, both experimentally and numeri- 
cally, to investigate Tes at various bath temperatures for 
various systems with varying degrees of anisotropy. 

In summary, I studied the fluctuation-dissipation ratio 
of the 3D Heisneberg SG with weak random anisotropy 
by off-equlibrium MC simulation. The calculated FDR 
exhibits a one-step-like behavior with the effective tem- 
perature of the SG state Tes — 2Tg irrespective of the 



bath temperature. The observed one-step-like behavior 
of the FDR is interpreted in terms of the recently devel- 
oped chirality scenario of SG transition. 

The author is thankful to M.Ocio, E.Vincent and 

H.Yoshino for useful discussion. The numerical calcu- 
lation was performed on the Hitachi SR8000 at the su- 
percomputer center, ISSP, University of Tokyo. 
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